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1. Introduction 

The calculation of the onium spectrum from first principles from lattice QCD represents a 
significant goal. There are a number of states that are stable to strong decay and far from thresholds 
(that might make finite size effects more significant). It should be possible to accurately calculate 
all of their masses. Further, since there are no valence up and down quarks, we have to consider 
only the sea quark mass dependence. However, one must be careful in dealing with heavy quarks 
on the lattice because aM is not small. Some early results using dynamical quark configurations 
were published in 2004 [|l|]. Using clover type quarks with the Fermilab interpretation we have 
been calculating the onium spectrum for some time ^ on MILC gauge configurations [^. The 
HPQCD/UKQCD collaborations have successfully been using NRQCD to treat the bottom quark 
on many of the same ensembles [Q]. Recently, they have started to use highly improved staggered 
quarks (HISQ) to study charm [|]. 

At Lattice 2005 |^], we presented results for lattice spacings a ?a 0. 18, 0. 12 and 0.09 fm. These 
are denoted extra coarse, coarse and fine, respectively. This year we have several new ensembles 
with a lattice spacing of 0.15 fm, denoted medium coarse in the plots. On these new ensembles, 
we have tuned the dynamical strange quark mass closer to its physical value based upon experience 
with the earlier ensembles. At the fine lattice spacing, we have new results on a more chiral ensem- 
ble with ami = 0.0031 and arus = 0.031 on a 40^ x 96 grid. All of our results for bottomonium are 
new. We also have some new results for the Xc2- (See Table [I] for details of the ensembles.) 

To calculate the onium spectrum, we need to find an appropriate value of the heavy quark 
hopping parameter K. We do this by studying the D., and Bs kinetic masses. We do this study on 
one ensemble for each lattice spacing and use the selected values of and K}, for all the ensembles 
with that lattice spacing. Errors in the kinetic masses tend to be large, and so we calculate mass 
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Table 1: Ensembles used in this calculation. 
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Figure 1: Splitting between the spin-averaged 
2S and \S states of charmonium. 
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Figure 2: Splitting between the Xdi'^P) and 
spin-averaged \S states. 



splittings based on differences in the rest energy. 



2. Charmonium Results 

We will examine a number of splittings in the charmonium spectrum. Let's look at the dif- 
ference between the spin averaged 2S and 15 states. In Fig. |l|, we plot the splitting in MeV as a 
function of the light bare sea quark mass. The experimental value is shown as a black fancy cross. 
The points in red are linear extrapolations in quark mass. We find that at each lattice spacing (ex- 
cept for the fine lattice) the extrapolated value is in good agreement with experiment. On the fine 
lattice, there is a considerable slope and the extrapolated value is quite high. (We note that as the 
horizonal axis is the bare quark mass, the slope is not a physical quantity; it also reflects a change 
in the mass renormalization as the lattice spacing changes.) 

Next we turn to issues of fine structure and look at the splitting between the ;^ci(l^) and the 
spin averaged 15 mass. It would be more natural to compare with the spin average of all the IP 
states, but the Xc2 mass is not yet available on all ensembles, so we compare with the IS spin 
average. Using the same color scheme and symbols as in the previous figure, we see in Fig. ^, 
that there is good agreement with experiment except on the fine ensembles which again exhibit 
a substantial slope. In this case, the two more chiral ensembles are in good agreement with the 
experimental value, but the ensemble with the heaviest value of the light sea quark mass gives too 
large a value for the splitting, leading to a large slope and too small a chiral limit. 

For the hc{\P) state shown in Fig. ^, we find good agreement with the experimental value on all 
ensembles. For the finest lattice spacing, there is only a modest slope from the chiral extrapolation. 

Turning next to the hyperfine splitting, we look at the 7/va(1S) - f]c{^S) mass splitting in 
Fig. ^. In this case, we find that the splittings are sytematically small, but that the value is increasing 
toward the experimental value as the lattice spacing decreases. 
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Figure 3: Splitting between the hc{lP) and spin- 
averaged 15 states. 



Figure 4: Hyperfine splitting of the 15 states. 
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Figure 5: Summary of charmonium spectrum. 



The Xcii^P) has only been studied on two ensembles so far. We have new results on one fine 
ensemble. In Fig. ^, we summarize the results for all the states studied. Except for the Xc2{iP), 
we plot results from our linear chiral extrapolation for each lattice spacing. For the ground states, 
if we focus our attention on the diamonds representing our smallest lattice spacing, we find the 
most serious discrepancy between our results and experiment is for the Xci- We have seen that 
our linear chiral extrapolation may be the culprit here, as the two more chiral ensembles are in 
good agreement with the experimental value. The S wave first excited states are not that well 
determined, but are rather heavy compared to the observed values. We have seen that on the finest 
lattice spacing, the high slope of the chiral extrapolation is accentuating the difference between our 
calculation and observations. Furthermore, the observed states are quite close to the DD threshold, 
which makes these states harder to calculate on the lattice without careful attention to finite volume 
effects. Thus, we are not seriously concerned about the high masses we are seeing for the 25 states. 
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Figure 6: Chiral extrapolation of the T(25)- 
T(15) splitting. 



Figure 7: Chiral extrapolation of the Xhi[^P) 
spin averaged 15 splitting. 



3. Bottomonium Results 

We have new results this year for bottomonium on the three smallest lattice spacings. Some 
important ground states in this system, the tj^ and hj,, have not been observed, so we shall modify 
some of the splittings we display. It is also interesting to compare our results with those using 
NRQCD for the heavy quarks [^. We expect that our results have larger discretization errors, 
when comparing them to the results of Ref. at a fixed lattice spacing, because the NRQCD 
action used in that work is more improved. 

In Fig. ^ we plot the splitting between the T(25) and T(15') masses. On the fine ensembles, 
the result is in good agreement with experiment. With larger lattice spacings, the splitting is a bit 
low. We look only at the T level since the T]/, masses are not well measured. In the bottomonium 
system, the 25 and 35 states are both below the BB threshold. Thus, possible finite size effects from 
a nearby threshold are not an issue, and we should get the 25 levels right. 

In considering the fine structure of the bottomonium spectrum we are faced with two issues. 
First, the hb has not been observed so we cannot compare our results with the experimental value 
(but we can make a prediction). Second, we have no results for Xbi yet, so we can't compute a spin 
average of the P states. So, we consider the Xhi, in the middle of the three ^Pj states. We compare 
it to the spin averaged 15 states assuming our r]b mass is correct. The result shown in Fig. 0has the 
splitting about 35 MeV too large on the fine lattices. 

The hyperfine splitting in the 15 states is shown in Fig. We find a value of about 38 MeV 
for this splitting for all a. The experimental result plotted comes from the PDG and is based on 
the observation of a single r]},. It should be noted that a preliminary result from CDF with higher 
statistics had a splitting of only 15 MeV. With this factor of 10 difference between the experimental 
results, it is hard to reach any definite conclusion about how accurate our hyperfine splitting is for 
the 15 states. It turns out, however, that on four ensembles we can directly compare our calculation 
with the HPQCD results. In Fig. ^, we compare the splitting in lattice units. The results are plotted 
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Figure 8: Hyperfine splitting of the bottomo- 
nium IS states. 



Figure 9: Comparison of HPQCD's hyperfine 
splitting using NRQCD with this work. 



on a semilogarithmic scale so that equal vertical differences represent equal fractional differences 
in the splittings. The upper two octagons and diamonds are coarse lattice results. The lower two 
octagons and diamonds are on fine ensembles. Points with the same bare quark mass are calculated 
on the same ensembles, but the HPQCD collaboration did not analyze every available configuration, 
so the actual set of configurations selected for analysis differs. We see that the fraction difference 
in the hyperfine sphttings is larger at the smaller lattice spacing and that the NRQCD hyperfine 
splitting is larger than that for clover. The difference is not unexpected, because the leading error 
for this splitting with the clover action is of order mv^, whereas in Ref. [^] it is of order mv^. 

In Fig. |l^, we summarize the spectrum of observed states using the spin-averged IS states 
(assuming a 38 MeV T-T]/, splitting) to set the additive constant needed to go from splittings to 
masses. Solid lines on the plot are experimentally determined masses. Dashed lines on the plot 
denote unobserved (or poorly observed) states. In the case of the Tj^, we use dashed lines and mark 
them CDF and PDG as discussed above. We find good agreement with experiment for the T(25'). 
The rih{2S) is in good agreement with theoretical expectations. Our masses for the P wave states 
hh, XhO and Xhi are too large. 

4. Conclusions and Outlook 

We observe only a mild dependence on the sea quark mass for most of the quantities studied 
here. In particular, the hyperfine splittings and the 25-15' splittings appear to be essentially inde- 
pendent of the sea quark mass. Hence a linear extrapolation seems adequate. For the P-wave states, 
we do observe a sea quark mass dependence at the 10% level on some ensembles. 

There are several positive features of the charmonium spectrum. The he and Xci masses look 
quite good (although the latter is driven low by the chiral extrapolation on the fine ensembles). 
At the smallest lattice spacing studied so far, the P wave splitting looks good (but see previous 
sentence). The hyperfine splitting is too small but improving as the lattice spacing decreases. 
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Figure 10: Summary of the bottomonium spectrum for ensembles with a k, 0.15, 0.12 and 0.09 fm. 

However, the 2S states are not accurately calculated. For the bottomonium spectrum, the excited 
state splitting looks good. It is not yet possible to test the hyperfine splitting, but our splitting is 
smaller than that coming from NRQCD calculations. Our P wave states seem too heavy. 

In the coming year, we expect to increase our statistics on a number of existing ensembles. In 
addition, MILC is generating ensembles with a lattice spacing of 0.06 fm and has plans to reduce 
the lattice spacing to 0.045. There is also some chance that there will be some ensembles with 
a = 0.105 fm. Currently, we are using an automatic criterion for picking the best fit. We need to 
consider alternative methods and whether picking a single fit properly reflects the systematic errors. 

The P wave splitting in the bottomonium spectrum does not seem to be in agreement with 
experiment. It will be interesting to see if the bottomonium spectrum improves as we reduce the 
lattice spacing (and am},). Kronfeld and Oktay [||] have been developing a highly improved clover 
quark action that we hope to use in the future. 
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